Observation of the Quantum Zeno Effect on a Single Solid State Spin 
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The quantum Zeno effect, i.e. the inhibition of coherent quantum dynamics by projective mea- 
surements is one of the most intriguing predictions of quantum mechanics. Here we experimentally 
demonstrate the quantum Zeno effect by inhibiting the microwave driven coherent spin dynamics 
between two ground state spin levels of the nitrogen vacancy center in diamond nano-crystals. Our 
experiments are supported by a detailed analysis of the population dynamics via a semi-classical 
model. 
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It was no later than 430 BC, when ancients greek's 
philosopher Zeno of Elea formulated his well known 
arrow paradox pQ. If an object can be observed at 
any instance of time to be at a well defined position in 
space motion can not occur. It took over 2 millennia 
until I. Newton and G.W. Leibniz solved the paradox 
with the invention of infinitesimal calculus. Then, in 
1977 it was shown that the paradox rises again from 
the formalism of quantum mechanics [2J. If a quantum 
system is permanently observed, its coherent dynamics 
can be significantly slowed down and eventually even 
frozen. Later on, the group of the 2012's Nobel laureate 
D.J. Wineland observed this so-called quantum Zeno 
effect in the radio frequency transition between two 9 Be + 
ground-state hyperfine levels [3]. Aside from philosophic 
implications of the phenomenon, recent proposals aim 
towards technologically utilizing the broader formulated 
quantum Zeno dynamics j4] for efficient quantum gates 
and fault tolerant quantum operation [5 7J. In this 
letter, we demonstrate the quantum Zeno effect in a 
solid state spin, namely the ground state spin of a single 
negatively charged nitrogen vacancy (NV) center in 
diamond. Beyond its analogy to a classical paradox the 
quantum Zeno effect here allows for a detailed study 
of the intricate interplay of coherent and incoherent 
dynamics of a single quantum system interacting with a 
macroscopic environment [5]. 

The NV center is one of the most studied color centers 
in diamond and single NV centers have proven to be ex- 
cellent solid state quantum systems. They are frequently 
used as single photon source [TO] , spin qubit [TflfTO] , 
or bio-compatible sensors [T3]. Nanodiamonds can be 
integrated into various other systems like cells, photonic 
crystal structures [ToT - flT] or plasmonic elements (181 119j. 
The origin of this versatility is a combination of the NV 
center's level structure, excellent ground state spin co- 
herence and stability of the diamond host lattice [TTl 120) . 
The electronic level structure according to the NV 
center's C3 V -symmetry is depicted in Fig. 1 (b) [2"T1 122j . 
The level structure consists of a triplet ground state 3 'A2 
with the spin m s = sub-level S Z A 1; and the m s = ±1 
states S x>y Ai. At zero magnetic field the m s = ±1 
states are split from the m s = level by approximately 
2.9 GHz. The triplet excited state 3 E has the two-fold 



degenerate m s = sublevels S z E Xiy , while the m s = ±1 
manifold is given by the four- fold degenerate S Xjy E Xjy . 
Using off-resonant excitation by a 532 nm laser levels 
of same spin can be linked by spin-preserving optical 
transitions at approximately 638 nm with the excitation 
rate k exc and radiative decay rate k ra d, assumed to be 
equal for all spin states. Furthermore, two singlet states 
1 A\ and 1 E exist. The state 1 Ai can be reached via inter 
system crossing (ISC) from the m s = ±1 manifold and 
decays via Ei to the m s = ground state as indicated 
by the dashed arrows in Fig. 1 (b). 

Via this ISC process the NV center is efficiently polar- 
ized within a few excitation-decay cycles. Furthermore, 
since the deshelving rate from Ex to the S Z A! ground 
state is comparably low the spin dependent ISC process 
makes the m s ± 1 states appear darker than the m s ± 
counterpart, enabling the optical detection of the 
magnetic resonance (ODMR). These spin detection and 
initialization mechanisms combined with long coherence 
times in the mainly spin-free diamond lattice render 
the NV center ideal to demonstrate coherent spin ma- 
nipulation by microwaves. For example electromagnetic 
induced transparency [23], simple quantum algorithms 
|24j and subdiffraction optical magnctometry |25j have 
been demonstrated. In Refs. 23. 25 a 532 nm laser is 




FIG. 1. (Color online) (a) The optical setup. Individual NV 
centers are excited by a pulsed 532 nm laser or a 532 nm 
cw laser modulated by an acousto-optic modulator (AOM) 
through a high NA objective lens. Fluorescence is separated 
by a dichroic mirror and an additional long pass filter (LP) 
prior to detection by a Hanbury-Brown and Twiss setup con- 
sisting of two avalanche photo diodes (APD) and a beam split- 
ter (BS). Additionally, microwave pulses can be applied via 
a thin gold wire in the vicinity of the diamond nano-crystals 
(MW) . (b) Schematic level diagram of the NV center without 
strain. See text for details. 



2 



used to inhibit a population transfer form the m s = 
to m s = +1 ground state, but it remains speculative, 
whether this is due to the quantum Zeno effect, or 
classical repumping into the m s = state. In this 
letter, we address this question experimentally, as well 
as by numerical simulations. For this purpose, at first 
the population dynamics of the investigated NV center 
are determined. In a second step, we use a pulsed 
532 nm laser to inhibit the coherent microwave driven 
population transfer between the m s = and m s = 1 
ground state. Our results are supported by numerical 
simulations based on the previously determined pa- 
rameters. Importantly, we demonstrate bidirectional 
inhibition of the population transfer, i.e. also for the 
transition from m s = 1 to m s = 0. In this case, 
classical repumping into the m s = state contra- 
dicts the effect, giving clear evidence for the existence 
of the quantum Zeno effect in this solid state spin system. 

In our experiment, we spin-coated commercial dia- 
mond nanocrystals with an average size of 50 nm (mi- 
crodiamant AG) onto a cleaned glass cover slip. Using a 
homebuilt confocal microscope (NA = 0.9) with Habury- 
Brown and Twiss setup a single NV center was identified. 
To perform the spin manipulation, microwaves pulses 
(MW) could be applied via a gold wire in the vicinity 
of the sample, while a small permanent magnetic field 
was used to split up the m s = ±1 spin levels by approx- 
imately 250 MHz. In the investigated nanodiamonds 
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FIG. 2. Fluorescence intensities Iq m (*) °f the NV center 
after switch-on of the cw laser when either the predominantly 
m s = state Go or the m s = 1 state Gi was initially prepared. 
The contrast defined by (7g — /gi)/(^g + ^Gi) is scaled up 
by a factor of 5 for better visibility. The solid curves are 
numerical solutions of the rate equations corresponding to 
the level diagram depicted in Fig. 3. 

magnetic splitting and transition rates vary for different 
centers due to strain and orientation. Thus, several mea- 
surements were performed to gain insight into the specific 
level-structure and dynamics. First, the well established 
ODMR technique [23J was used to identify the ground 
state spin resonances. Subsequently, coherent Rabi oscil- 




FIG. 3. Simplified level diagram and transition rates of the 
NV center with strain. The m a = (1) electronic ground 
state is denoted Go m , while the corresponding predominantly 
m s = 0(1) excited states are E (i). The singlet states are 
merged to the state S. Solid arrows correspond to allowed 
transitions, while dashed arrows indicate forbidden transi- 
tions. The microwave induced coherent ground state spin 
rotations are indicated by MW. The inset shows the pulse 
sequence used to measure the Zeno effect. MW denotes the 
microwave source used to apply up to 3 7r-pulses. The 532 nm 
cw laser is used to initialize the NV center at the beginning 
of the experiment and for state readout at the end, while a 
strong pulse from the pulsed laser is applied at time delay r 
with respect to the center of the central 7r-pulse. The NV 
fluorescence intensity is measured within short time windows 
after and prior to applying the 7r-pulses to obtain the signal 
and normalization reference as indicated. 



lations were driven to measure the oscillation frequency 
f2, as well as the damping of the oscillation I/T2 (cf. 
Supplemental Material) . To estimate the transition rates 
between different levels, first the m s = state was pre- 
pared by applying a green cw laser (1.2 mW) for of about 
5 /xs (cf. Fig. 1 (a)). About 2 /is after switching off the 
laser the NV center is assumed to be relaxed into the 
desired m s = state, which can be transformed into the 
m s = 1 state by an optional MW 7r-pulse. Subsequently 
the time dependent fluorescence intensity after switching 
on the cw laser again was measured (Fig. 2). These dy- 
namics strongly depend on the transition rates between 
different levels. To estimate these transition rates, we 
have to remark that strain induces a mixing of exited 
states with different spin orientations by the angle 1?, re- 
sulting in the model depicted in Fig. 3, where we denote 
the predominantly m s = 0(1) excited state Eq(i). The 
mixing has two direct implications: Spin non-preserving 
optical transitions, as well as inter-system crossing from 
the predominantly m s = excited state become possi- 
ble [2S1I2Z]- According to Fermi's golden rule the rates 
of the transition allowed in the unstrained NV are re- 
duced by the factor cos($) 2 , whereas the formally for- 
bidden transitions have now the rates k x sin(i9) 2 , with 
k x being the rate of the corresponding allowed transi- 
tion. To simplify the model, we do not consider the 
m s = — 1 manifold and merge the singlet states to the 
single deshelving state S. A differential equation repre- 
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sentation of the model can be found in the Supplemental 
Information. Using this model the measured dynamics 
can be exactly reproduced (solid line in Fig. 2). A sum- 
mary of the deduced rates is given in Table 1. A first 
important consequence of the model is a reduced polar- 
ization efficiency i] po i. By illumination with green light 
the bright state G = r] po iG + (1 - r) po i) exp(iy>i)Gi 
with random phase ipi is prepared. A subsequent mi- 
crowave 7r-pulse transfers the population into the dark 
state Gi = 77 po /Gi + (1 — rj po i) exp(«(^2)Go- For additional 
verification of the deduced parameters, we independently 
measured the excited states lifetimes with a pulsed laser 
(PicoQuant LDH-P-FA-530-XL), deduced the excitation 
rate k exc from saturation measurements and the derived 
contribution of fluorescent background to the signal h g 
from autocorrelation measurements. All parameters were 
consistent with those deduced from the fit. 



TABLE I. Transition rates and parameters deduced by fitting 
the model (Fig. 3) to the measurement shown in Fig. 2. 
The error corresponds to one standard deviation confidence 
interval. 





(80.0 ± 0.5) ns 


T 2 


(1.0 ±0.5) MS 


1 / kexc 


(33.5 ± 5)ns 




(34±5)ns 




(146±9)ns 


l/k s 


(38±5)ns 




(20.1 ±1.7)° 


Ibg 


0.09 ±0.09 


rjpol 


0.85 ±0.01 



After determining the transition rates the quantum 
Zeno experiment was simulated and experimentally re- 
alized. The experimental sequence is illustrated in the 
inset of Fig. 3. First, the NV center is initialized to the 
bright state Go by applying the green cw laser for about 
5 /is. In order to initialize the dark Gi state a subse- 
quent MW pulse can be applied for a length of 40 ns. 
This corresponds to a 7r-pulse and transfers the popu- 
lation from Go to &x. After the initialization the MW 
pulse is switched on starting a coherent transition from 
Go to Gi (or from Gi to Go). The pulse is set to a fixed 
length of 40 ns, i.e. a 7r-pulse. Synchronized to the mi- 
crowave pulse, a strong green laser pulse from the pulsed 
laser (pulse energy about 350 pJ) is applied at varying 
time delay r. This laser pulse performs a projective mea- 
surement of the NV center's spin state, although the gen- 
erated photons are not detected. Finally, about 300 ns 
after the microwave pulse the green cw laser is turned on 
again and the fluorescence of the NV is recorded, giving 
a measure of the remaining population in the brighter Go 
state. Fig. 4 shows the detected intensity (dots) and nu- 
merical simulations (solid). Dashed vertical lines indicate 



the switch-on and switch-off of the MW pulse. Quantum 
theory predicts, that the strong projective laser pulse de- 
stroys the microwave induced coherent polarization, and 
thereby inhibits the dynamics. This process depends on 
the time delay between the laser and the MW pulse. At 
the center of the MW-pulse the polarization reaches its 
maximum (cf. Supplemental Material). Hence a pro- 
jective laser pulse at r = effectively inhibits further 
coherent dynamics and the final state has a large compo- 
nent of the initial state. This behavior is clearly visible 
in Fig. 4(a) as increased (decreased) fluorescence around 
r = when the initial state was Go (Gi). The experi- 
mental results are well reproduced by the numerical sim- 
ulation based on the previously determined parameters. 
In particular, the decreasing intensity for an initial Go at 
positive values of r in Fig. 4 (a), i.e. when the NV is in 
the dark state Gi when the strong laser pulse is applied, 
proves that the effect is not due to repumping of the NV 
center. 

read-out of G-, 




-50 50 
x in ns 



-50 50 
x in ns 



FIG. 4. Measurement (dots) and numerical simulation (solid 
lines) of the NV fluorescence after the quantum Zeno experi- 
ment, (a) The NV center is initialized to the bright Go (dark 
Gi ) state. The coherent population transfer during the MW 
pulse (indicated by dashed lines) is inhibited by a strong green 
laser pulse at time delay r with respect to the center of the 
MW pulse. Subsequently, the fluorescence as a measure of the 
occupation of the bright state Go is probed. As predicted by 
the numerical simulation (solid curve), the Zeno pulse is most 
effective at r = 0. (b) Same as (a), but with an additional tt- 
pulse before measuring the NV fluorescence, i.e. monitoring 
of the population of the dark state Gi . 

To further support this interpretation, we repeated 
the experiment with an additional 7r-pulse before mea- 
suring the NV fluorescence. This swaps the population 
of the dark and bright state, i.e. the measured inten- 
sity is proportional to the population of the dark state 
Gi (Fig. 4 (b)). Again, as predicted by the simulation 
(solid curve), the transition is most effectively inhibited 
when the projective laser pulse is applied in the center 
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of the microwave 7r-pulse. Here the curve for an initial 
Gi, where the coherent spin rotation from Gi to Go is 
inhibited by the laser, is particularly remarkable. In this 
configuration the laser pulse effectively increases the pop- 
ulation of the Gi state, resulting in an increased fluores- 
cence intensity when probing. To get further insight, we 
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FIG. 5. Experimental values (dots) and numerical simula- 
tion (solid curves) of the m s spin projection when initially 
preparing the predominantly m s = (1) level. Time inter- 
val between the vertical lines corresponds to the MW pulse, 
while the shaded areas are inaccessible due to imperfect po- 
larization. The features centered at r = are whiteness of 
the Zeno effect. Dashed lines indicate the expectations for 
the ideal case, where the laser pulse destroys the coherence 
without altering the populations. 

calculated the spin projection < m s > from the contrast 
between the intensity without final 7r-pulse (Fig. 4(a)) 
and the intensity with final 7r-pulse (Fig. 4(b)). After 
preparation of the initial state the polarization efficiency 
is r/poi = 0.85 as found from fitting the data (Fig. 2) 
to the model (Fig. 3). Therefore we set < m s >= 0.15 
to correspond to the contrast after preparation of Go, 
while < m s >= 0.85 corresponds to the contrast after 
preparation of Gi. For the ideal case, where the laser 
pulse destroys the coherence without altering the popu- 
lations, a laser pulse not overlapping with the MW pulse 
should have no influence, i.e. the initial population is 
completely swapped. A laser pulse within the MW pulse 
inhibits the coherent dynamics so that at r = at most 
50% population swapping is reached. The deviation of 
measurement and realistic simulation form this ideal case 
(cf. Fig. 5) has several reasons. These are a limited Ti 
time, repolarization, and a finite lifetime of the excited 
states. The first one damps the Rabi oscillations, re- 
sulting in a small offset towards < m s >= 0.5, partic- 
ularly visible at the curve for an initial Go for negative 
r and for positive r at the curve for an initial Gi. Re- 
polarization is effective when the NV state has a large 
Gi component at time r. This drives the spin towards 
< m s >= for positive (negative) r in case of initial 
preparation of Go (Gi). The finite lifetime of the NV ex- 
citation is an issue, when the laser pulse is applied before 



the end of the microwave pulse, i.e. for r < 20 ns. The 
laser pulse drives the NV to the excited states, where 
the MW pulse is ineffective. Consequently this effect is 
most relevant when the laser pulse is applied immedi- 
ately before the MW pulse at r = —20 ns, while it gets 
less relevant for smaller or larger r. The behavior is gov- 
erned by the excitation lifetime r exc and is proportional 
to 1 — cxp(t/t exc ) for r < —20 ns. Within the MW pulse 
(—20ns < t < 20ns) the shape is given by exp(r/r ea;c ). 
For initial Go the excitation lifetime r exc is on the order 
of r exc sa l/k ra d while for an initial Gi the branching ra- 
tio into the metastable singlet S is much larger giving a 
excitation lifetime T exc « 2/(k ra d + kd es h)- This influence 
of the finite excitation lifetime can well explain the de- 
crease at r < —20 ns in the curve for initial G , but does 
neither explain the shape nor position of the feature at 
t = 0. This feature is caused by the destruction of the 
quantum coherence by the strong measurement pulse, i.e. 
the Zeno effect. 

In conclusion we have demonstrate that a strong laser 
pulse represents a projective measurement of the spin- 
state of single NV center. A microwave-driven coherent 
time evolution can be inhibited in this way, representing 
the quantum Zeno effect in a solid-state system. A model 
including realistic and experimentally confirmed param- 
eters was derived to explain the quantum Zeno dynamics 
in detail. Future work will be devoted to the demon- 
stration of the quantum Zeno effect under spin selective 
optical excitation, e.g. at cryogenic temperature. Here 
the effect will not only be more pronounced, but is also 
relevant for recently proposed robust two-quibit quantum 
gates 
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SUPPLEMENTARY INFORMATION 




FIG. 6. Measured Rabi oscillations. A resonant microwave 
pulse of length t is applied and the fluorescence intensity 
recorded. The solid line is a fit to the data. MW power was 
slightly adjusted for the Zeno experiment to reach Q — 80 ns. 
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FIG. 7. Calculated Rabi oscillations of an ideal two level 
system without damping. Top: Calculated Go(t) and Gi(t) 
occupation during a 3n microwave pulse when starting with 
Go(0) = I. Bottom: The corresponding microscopic polariza- 
tion Goi(t) reaches its maximal values at Qt = (n + 1/2)tt. 
Then, a projective measurement most efficiently inhibits the 
dynamics. 



The time derivatives of the level populations (cf. main text Fig. 3) and the quantum mechanic polarization between 
the ground states Goi(i) are given by Eqn. (1-6). Here, il(t) is the time dependent MW Rabi frequency, which is equal 
to n during the MW pulses and zero else. The excitation rate is modulated in a similar manner. Here we assume the 
switch-on to be of Gaussian shape with a HWHM of 2.5 ns, which is consistent with the measured slope of the used 
AOM. The pulsed laser is simulated by a 1 ns long interval of k exc = 2.5/ns. The equations were numerically solved 
using Wolfram Mathematica, while initial values Gq(0) and Gi(0), i.e. Go and Gi were obtained to be self consistent. 



E (t) = k exc (t) [aGo(i) + /8Gi(t)] - [k rad + pk s ] E (t) 
Exit) = keccdt) [oGi(t) + (3G (t)} - [k rad + ak s ] E±{t) 

1 



Goi(t) = --m [G x {t) - G (t)] - Goi(t) 



+ k exc (t) 



G (t) = -Im[n(t)Goi(t)] - k exc (t)G [t] + k rad [aE {t) + pE x {t)] 
G x {t) = Im [fi(t)G i(t)] - fc exc (t)Gi[t] + k rad [aExit) + pE (t)} 
S(t) = k s \aE x {t) + pE {t)\ - k desh S(t). 



k des hS(t) 



(1) 
(2) 

(3) 

(4) 
(5) 
(6) 



r 



